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ABSTRACT 

We present multiple Chandra and XMM-Newton observations of the type 1.8 
Seyfert Galaxy NGG 1365, which shows the most dramatic X-ray spectral changes 
observed so far in an AGN: the source switched from reflection dominated to 
transmission dominated and back in just 6 weeks. During this time the soft 
thermal component, arising from a ~ 1 kpc region around the center, remained 
constant. The reflection component is constant at all timescales, and its high flux 
relative to the primary component implies the presence of thick gas covering a 
large fraction of the solid angle. The presence of this gas, and the fast variability 
time scale, suggest that the Gompton-thick to Gompton thin change is due to 
variation in the line-of-sight absorber, rather than to extreme intrinsic emission 
variability. We discuss a structure of the circumnuclear absorber/reflector which 
can explain the observed X-ray spectral and temporal properties. 

Subject headings: Galaxies: AGN — X-rays: galaxies — Galaxies: individual 
(NGG 1365) 


1. Introduction 

According to the Unihed Model of Active Galactic Nuclei (AGNs, see review by Urry & 
Padovani, 1995), an axisymmetric absorber/reflector is present around the central black hole 
of AGNs, with a size between that of the Broad Emission Line Region (BELR, ~ 1000 Rs, 
Schwarzschild radii, corresponding to ~ 10“^ pc for a 10® Mq black hole) and that of the 
Narrow Emission Line region (tens to hundreds parsecs). The simplest geometrical and 
physical conhguration of such absorber is that of a homogeneous torus (Krolik & Begelman 
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1988) beyond the dust sublimation radius {vgub ~ Barvanis 1987, where Ljjv,m is 

the UV luminosity in units of 10^® erg s“^). However, this view has been recently challenged 
by several pieces of observational evidence: 

(1) Dramatic X-ray absorbing column density changes (factors of >10) over a few years have 
been seen in several type 2 (narrow permitted line) Seyferts (Risaliti, Elvis & Nicastro 2002), 
ruling out a homogeneous absorber. 

(2) Rapid column density variability, on time scales of hours, requires an X-ray absorber no 
larger than the BELR. Such changes have been detected in the brightest absorbed Seyfert 
Galaxy, NGC 4151 (~ 10 — 30 hours, Puccetti et ah 2004) and in the Seyfert 2 Galaxy 
NGG 4388 (4 hours, Elvis et ah 2004). Assuming that the absorber is made of material 
moving with Keplerian velocity around the central source, these observations imply that its 
distance is of the order of that of the broad line clouds, i.e. ~ 10^ i? 5 , where Rs = 2 GMbh / 
is the Schwartzschild radius. 

(3) The reflection components in the X-ray spectra of Gompton thin Seyfert Galaxies are 
systematically stronger than expected from reflection off gas with the same column density 
measured in absorption. This has been shown convincingly for a few sources in which a 
detailed measurement of the reflection component has been possible (NGG 1365, Risaliti et 
ah 2000; NGG 2992, Gilli et ah 2002; NGG 6300, Guainazzi et al. 2003) and in a statistical 
sense in a sample of ~ 20 bright Seyfert 2s (Risaliti 2002). 

NGG 1365 (;2=0.0055) is a particularly striking example of extreme X-ray variability 
and highly efficient reflection. It was observed by ASGA in 1995 (lyomoto et al. 1997) 
in a reflection-dominated state (Thomson optical depth r > 1, corresponding to N^^ > 
1.5 X 10^^ cm“^) and, three years later, in a Gompton-thin state {Nh ~ 4 x 10^^ cm“^) by 
BeppoSAX (Risaliti et al. 2000). The 2-10 keV flux of the reflection emission was more than 
5% that of the intrinsic spectrum measured by BeppoSAX. This corresponds to a ratio R=1 
between the normalizations of the reflected and direct components in the PEXRAV reflection 
model (Magdziarz & Zdziarski 1995) Such a high reflection efficiency can be achieved only if 
a thick (Nh > 3 x 10^^ cm“^) reflector covers a large fraction of the solid angle around the 
central source (Ghisellini, Haardt & Matt 1994). 

Similar switching from reflection-dominated to transmission dominated states has now 
been observed in a handful sources on timescales of a few years (Matt, Guainazzi & Maiolino 
2004) and could be due either to extreme variations of the intrinsic luminosity and so of the 
ionization state of the absorber (as suggested by Matt et al. 2004), or to Gompton-thick 
clouds moving across the line of sight, so changing the line of sight column density^. 


^It is worth noticing the conceptual difference between “reffection dominated” (which merely refers to 
the observed spectral shape) and “Compton thick” (which implies a physical interpretation of the observed 
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Here we present an X-ray observational campaign on NGC 1365, consisting of 3 XMM- 
Newton observations, plus a Chandra observation (Obsid 3554) performed just three weeks 
before the first XMM-Newton observation. The observing times are relatively short (~ 
10 — 20 ksec), therefore a temporal analysis within single observations is not possible. We 
will therefore discuss the spectra extracted from the whole observations. Two more, longer 
(60 ksec), XMM-Newton observations were recently performed, as a continuation of our 
program of monitoring of this source. Detailed spectral analysis of all the XMM-Newton 
observations will be presented in a forthcoming paper (Risaliti et al. 2005, in prep.) 

Here we will concentrate on the analysis of the most striking result of our study, i.e. dra¬ 
matic hard X-ray continuum variations: changes from reflection-dominated to transmission- 
dominated spectra occurred on times scale shorter than three weeks. During this time, the 
soft thermal component remained constant in all observations. The high resolution of Chan¬ 
dra allowed us to resolve this emitting region, which extends over ~ 1 kpc from the center, 
while the hard component originates in a region < 200 pc (2 arcsec) diameter^. 


The observation log of the X-ray observations of NGC 1365 is reported in Table 1 
together with basic spectral htting results (§2). 


2. Data Analysis and Results 

The spectra presented here have been obtained with CCD detectors: the ACIS-S in¬ 
strument on board Chandra, and the EPIC PN and MOS instruments on board XMM- 
Newton. All observations were performed in full-frame mode. We did not find any signif¬ 
icant pile-up. This was expected for XMM-Newton (the previous published 2-10 keV flux 
is < 10^^ erg s“^ cm“^, mostly emitted at energies > 3 keV) , while for Chandra this is a 
consequence of the extremely faint state of the source (Table 1) at that time. 

The data were reduced using the standard procedures, using the CIAO v.3.0 package 
for the Chandra observation and SAS v6.0 for the XMM-Newton observations. The spectra 
were extracted from circular regions of a radius of 2” ( Chandra) and 30” (both Chandra and 
XMM-Newton). The background spectra were extracted from regions close to our source 
and free from contamination by bright serendipitous sources. The background contribution 
is negligible for the Chandra spectrum of the central 2” radius region (thanks to the small 


reflection dominated spectrum as due to Compton thick absorption of the direct component). 


^Throughout this paper we adopt Hq = 70 km s ^ Mpc ^ (Spergel et al. 2003). 
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extraction region), and for the two XMM-Newton spectra when the source was caught in a 
bright state (XMM 1 and XMM 3). In the remaining two spectra (XMM 2, and the Chandra 
spectrum from the 30” radius region) the background contamination is ~ 10%. 

The spectra were analyzed using the XSPEC 12.0 package. The spectra of the XMM- 
Newton observations are shown in Fig. 1. A visual inspection of this hgure is sufficient to 
clearly notice the main spectral variation: in the 2-10 keV band the XMM 2 spectrum is an 
order of magnitude weaker than the XMM 1 and XMM 3 observations and has a prominent 
iron line with an equivalent width EW> 1 keV. The continuum is also signihcantly flatter 
than in the XMM 1 and XMM 3 5-10 keV spectra (at lower energies, the XMM 1 and XMM 3 
spectra are inverted, due to a photoelectric cut off). These features strongly suggest that 
during the XMM 2 observation the source was in a reflection-dominated state. 

One more interesting aspect is the constancy (within 5%) of the soft emission. This 
suggests that the soft spectrum is dominated by a component not directly related to the 
primary AGN emission. This is conhrmed by the visual analysis of the Chandra spectra 
of the central 2” region, and of the 30” region; most of the soft emission comes from the 
extended region between 2” and 30” (corresponding to a sphere of outer radius of ~ 2 kpc). 

In order to conhrm the results of this visual inspection, we performed a detailed spectral 
analysis of each observation. The complete report on this work will be presented elsewhere 
(Risaliti et ah 2005, in prep.). In the following we briefly summarize the results relevant for 
our interpretation of the continuum variations. 

Chandra. The G/iandra spectrum obtained from the central 2” is unusually flat: the photon 
index of a simple power law model is T = —0.1 ± 0.1. A good £t is obtained with a cold 
reflection model (PEXRAV in XSPEC), plus an iron emission line at 6.4 keV (equivalent width 
EW=200+}^° eV). A soft thermal component is not required by the data. The spectrum 
obtained from the larger 30” radius region can be fitted with the same components, plus 
thermal emission with kT~0.8 keV, with a normalization consistent with that found in all the 
XMM-Newton observations. This implies that a) the hard emission comes from the central 
region; b) the soft component originates from a larger nuclear region (radius ~ 1 kpc). 

XMM 1. The first XMM spectrum, obtained only 3 weeks after the Chandra observation, 
is clearly dominated by the direct emission of the AGN above ~ 3 keV. The best fit model 
consists of all the components used to fit the Chandra spectrum, plus an absorbed power law, 
completely dominating the emission above a photoelectric cut off at 3keV{NH = 4.8;[|;^x 
10^^ cm“^) The best fit parameters for the cold reflection are fully consistent with those 
found in the Chandra observation. The ratio between the normalizations of the reflection 
component and the power law is R=1.2. 
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XMM 2. The second XMM spectrum was obtained 3 weeks after the hrst. The spectrum is 
now back in a reflection-dominated state (Fig. 1). The best £t model consists of two compo¬ 
nents: a cold reflection emission, and an iron emission line at E=6.4 keV, with EW~1.2 keV. 
All the parameters of the cold reflection component are consistent with those found in the 
previous XMM 1 and Chandra spectra. 

XMM 3. This spectrum, obtained 5 months after the XMM 2 observation, is dominated by 
the direct emission from the central source. The best £t model is the same as in XMM 1, with 
a normalization of the direct component higher by 50%, and an absorbing column density 
Nh = 3.41q2 X 10^^ cm“^. The reflection component is consistent with that of XMM 1 and 
XMM 2. The ratio between the normalizations of the reflection component and the power 
law is R=0.8. 

The fit is significantly improved (Ay^ = 45) by the addition of a broad (possibly rela¬ 
tivistic) iron emission line^, although other models, e.g. partial covering and a second, lower 
Nh, absorber (as in NGC 4151, Zdziarski et al. 2002) may also fit. 


3. Discussion 

The spectral analysis of multiple Chandra and XMM-Newton observations of NGC 1365 
reveals transitions from reflection-dominated to transmission-dominated, on time scales no 
longer than 3 weeks. This is, by a factor of 30, the shortest time scale observed so far for 
such extreme X-ray spectral variations (Matt et al. 2003) While we also find a kiloparsec 
scale diffuse origin for the soft thermal component and the possible presence of a relativistic 
iron line in the highest statistics observation, we will defer discussing of these issues to where 
we can also consider two new, longer (~ 60 ksec) XMM-Newton observations (Risaliti et al. 
2005, in prep.). Here we concentrate on the fast spectral transitions. 

The extreme column density variations observed from the Chandra and XMM 1 observations, 
and from XMM 1 and XMM 2 can be due to either: (A) fading of the central source, down 
to a flux 50 times or more fainter than in the “transmission dominated” observations; (B) 
increased ionization making the absorber more transparent; (C) column density variability, 
with the reflection dominated spectra having Nh > 10^^ cm“^ along the line of sight, 

(A) Fading of the central source. In the intrinsic variability scenario, an almost 
complete switch-off of the source (down to a flux < 2% of the intrinsic flux in the XMM 1 


^This is an important point in itself, however it does not affect our treatment of the continuum variations, 
therefore we will not discuss it further here, but in Risaliti et al. (2005). 
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observation) is needed, in a time shorter than 3 weeks. This can be compared with the 
typical cooling time of the inner part of the accretion disk^. The thermal timescale in a 
Shakura-Sunyaev (1972) disk is approximately tth ~ (a x (Frank, King & Raine 

2002, Starling et al. 2004), where a is the viscosity parameter [a < 0.1 for a realistic 
model), and fix is the Keplerian angular velocity, fix = ^/GMbh/R^- The black hole mass, 
Mbh can be estimated from the K magnitude of the galaxy (K=8.4, Jarrett et al. 2003), 
through the relation between black hole mass and luminosity in the K band (Marconi & 
Hunt 2003). We obtain Mbh ~ 1.5 x 10® M©, and tth > 30 i? 2 o^ days, where R 20 is the 
linear dimension of the X-ray source in units of 20 Rs- Given the uncertain assumptions this 
result is in marginal agreement with our observational requirement t < 21 days. However, 
this requires an extreme scenario, i.e. the complete switch off of the source right after the 
XMM 1 observation. We conclude that intrinsic variability due to disk cooling is an unlikely 
explanation for the observed spectral changes. However, in order to rule it out completely, 
shorter variations (in a timescale of a few days) should be observed. 

(B) Ionization changes. In order to make the cold absorber found in the Chandra and 
XMM 2 observations transparent, an increase of the ionization parameter by at least ~ 3 
orders of magnitude would be needed (from .^ < 1, as required in order to completely obscure 
the 1-10 keV emission, to .^ > 10®, as required in order to make the gas transparent. This 
would in turn imply a similar increase in the intrinsic luminosity, which is not seen. 

(C) Nb variability. We assume (as in Risaliti et al. 2002, Elvis et al. 2004) that the 
absorber is made up of gas clouds moving around the central source with Keplerian velocity 
(motions close to Keplerian dominate in the high ionization parts of the BELR in AGNs, 
Peterson & Wandel 2000). The variability timescale is then given, to a hrst approximation, 
by the crossing time for a cloud across the line of sight. Assuming spherical clouds, the 
distance R of the absorber from the center is (Risaliti et al. 2002): R 600^100^10-^^,24-^5 
where tioo is the variability time in units of 100 ksec, nio is the cloud density in units of 
10^° cm“® and Nh, 2 a is the column density of a single cloud in units of 10^^ cm“^. Assuming 
that the variations from Gompton thick to Gompton thin states were due to a single cloud 
passing along the line of sight, the measured Nh variation, A(A’x) = 10^^ cm“^, corresponds 
to the column density of a single cloud. We further assume the cloud crossing time to be 
no longer than the time interval between the XMM 1 and XMM 2 observations, fioo A 20. 
With these numbers we obtain R < 2 x lO^nf^Rs- A further stringent requirement for the 
absorbing cloud is to be large enough to cover the X-ray source. Using Mbh = 1.5 x 10® Mq 


"‘Here we assume that the fading in the X-rays is due to a decrease of the seed photons from the disk 
which are up-scattered by the coronal hot electrons. 


= L{2 — 10 keV)/nR^, where n is the number density of the absorbing gas. 
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(see above), and assuming again a linear dimension of the X-ray source D~20 Rs, we obtain 
n < 2 X 10® cm“^ and R < 10^ Rs = 4.5 x 10^^ cm. For comparison, the dust sublimation 
radius (assuming the standard X-ray to bolometric correction of Elvis et ah 1994) is of the 
order of 10^^ cm. 

The Compton thick - Compton thin variation discussed here is, by a factor ~30, the fastest 
observed so far (Matt et al. 2003, Guainazzi et al. 2004). Therefore, the arguments used to 
rule out the intrinsic variability scenario are significantly less stringent in the other known 
cases. As a consequence, it is at present impossible to tell whether the observed extreme 
variability in NGC 1365 is an unique case or is representative of the other known state¬ 
changing sources. However, several general considerations can be made on the occurence of 
such variations, which are relevant in modeling the circumnuclear medium of AGN, and in 
driving further studies: 

- Compton thick gas seems to be present in (almost) all X-ray obscured Seyfert Galaxies. 
In ~40-50% of the sources it covers the line of sight to the central X-ray source (Risaliti, 
Maiolino & Salvati 1999). In the other X-ray obscured Seyfert Galaxies bright enough to 
allow a careful X-ray spectral analysis, the high reflection efficiency suggests the presence of 
Compton thick gas reflecting the primary emission (see Introduction). 

- Even if we ascribe all the known thick-thin transitions to Nh variability, this event is still 
quite rare: only four objects are known to show such variations (Matt et al. 2003) out of 
the more than 30 known Compton thick AGNs (Comastri 2004). Recently Guainazzi et al. 
(2004) found one more example (NGC 4939) in a sample of 11 objects. 

- The few indications available on short time scale variations (for example this work, and 
Elvis et al. 2004) suggest that the absorber is rather compact, at a distance of 10^ — 10^ Rs 
from the central source. 

This observational evidence can be explained if we assume a stratihed structure for 
the circumnuclear absorber/reflector, with a central planar Compton-thick region and a 
decreasing average column density at greater angles or distance above the disk (Fig. 2). The 
gas is clumpy and at the distance of the order of that found in NGC 1365 (~ 10^ Rs). 
The average number of clouds along a given line of sight is A^ ~ 5 — 10 (except, possibly, 
for the completely Compton thick region), in order to reproduce the observed average 
variability in bright AGNs (Risaliti et al. 2002). The angles have been chosen in order to 
reproduce the observed column density distribution of Risaliti et al. (1999), with fractions of 
45%, 25% and 10% of the solid angle covered by gas with column density Nh > 10^^ cm“^, 
Nh ~ 10^^ —10^^ cm“^, and Nh ~ 10^^ —10^^ cm“^, respectively, and with the remaining 20% 
free from absorption, in agreement with the obscured/unobscured Seyferts ratio of Maiolino 
& Rieke (1995). 

In this scheme strong absorption variations are unlikely in the extreme cases of heavily 



Compton-thick sources (line of sight through the central region of the absorber) and low 
obscuration [Nh = 10^^ — 10^^ cm“^) sources (line of sight far from the absorber), while 
they can happen if the line of sight intersects the transition region between the Compton 
thick and Compton thin zones. In this simple scheme, NGC 1365 is one of the sources seen 
through this “transition region”. We note that a similar stratihed and compact structure 
has been proposed to explain the relative X-ray/mid IR absorption properties of Luminous 
Infrared Galaxies (Risaliti et al. 2000), and to explain the Nh distribution and the megamaser 
emission in AGNs (Kartje, Kbnigl & Elitzur 1999). In this latter case the distribution arises 
as the result of a disk wind. 

The proposed scheme needs to be tested with more detailed studies of the single sources 
showing the most extreme variability (a new XMM-Newton observational campaign is on 
going) and with an unbiased statistical analysis of the occurence of such extreme variations. 


We are grateful to A. Siemiginowska and G. Matt for useful discussions. This work was 
partially supported by NASA grants NAG5-13161, NNG04GF97G, and NAG5-16932. 
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Fig. 1.— Upper panel: Unfolded spectra (nsing the best fit models of Tab. 1) of the XMM 1, 
XMM 2 and XMM 3 observations. Central Panel: nnfolded spectra of the Chandra obser¬ 
vation, extracted from a circnlar region of 2 arcsec (circles) and of 30 arcsec (triangles), as 
for the XMM observations. Fluxes are in units of 10“^ erg s“^ cm“^ keV“^. 
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Fig. 2.— Section of the proposed structure of the circumnuclear absorber/reflector. The 
average column density is indicated by the colors of the background and of the clouds (dark: 
Nh > 10^^ cm“^; medium: 10^^ cm“^ < Nh < 10^"^ cm“^; light: 10^^ cm“^ < Nh < 
10^^ cm“^). The angles are chosen in order to reproduce the column density distribution of 
Risaliti et ah (1999). 








Table 1: NGC 1365 - X-ray Observation log. & Spectral Fitting Results 


Observatory 

Date 

Counts 


F 



F(2-10)‘^ 

1 

o 

xl 

Chandra'^ 

2002 Dec 24 

7478 

14.6 

i-98^S:S(g) 

- 

- 

0.21 

0.14 

1.09 

XMM-1 

2003 Jan 16 

14866 

17.8 



1.2 

0.48 

1.23 

1.09 

XMM-2 

2003 Feb 09 

3503 

5.8 

2.09«:!l(g) 

- 

- 

0.17 

0.11 

1.06 

XMM-3 

2003 Aug 13 

9466 

8.7 

o 00+0.12 

33.5l^;9 

0.8 

0.70 

2.30 

1.21 


Duration (ksec). in units of 10^^ cni“^. Ratio between the normalizations 
of the reflected and transmitted components. Observed 2-10 keV flux, in units of 
10“^^ erg s“^ cm“^. 2-10 keV luminosity in units of 10"^^ erg s“^. For the transmis¬ 

sion dominated spectra, the luminosity is absorption corrected. 2 arcsec dia. aperture. 
photon index of the intrinsic power law in the PEXRAV model. 



